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Abstract. 

It is shown that the hydrodynamic modes of a dilute granular gas of inelastic hard 
spheres can be identified, and calculated in the long wavelength limit. Assuming they 
dominate at long times, formal expressions for the Navier-Stokes transport coefficients 
are derived. They can be expressed in a form that generalizes the Green-Kubo relations 
for molecular systems, and it is shown that they can also be evaluated by means of 
A-particle simulation methods. The form of the hydrodynamic modes to zeroth order 
in the gradients is used to detect the presence of inherent velocity correlations in the 
homogeneous cooling state, even in the low density limit. They manifest themselves in 
the fluctuations of the total energy of the system. The theoretical predictions are shown 
to be in agreement with molecular dynamics simulations. Relevant related questions 
deserving further attention are pointed out. 
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1. Introduction 

Rapid granular flows exhibit a phenomenological behavior that is similar in many ways 
to that of ordinary fluids. Then it is not surprising that the methods of statistical 
mechanics and kinetic theory are being successfully extended to describe them. The 
primary feature of grains that must be captured in any idealized model is the inelasticity 
of collisions. The prototype is a system of smooth inelastic hard spheres (or disks in two 
dimensions) whose energy loss in collisions is characterized by a coefficient of normal 
restitution. In real granular media, this coefficient depends on the relative velocity of 
the colliding particles pQ , but for the sake of simplicity, and expecting to retain the main 
qualitative features, it will be considered here as a material constant. 

At the most fundamental level of particle dynamics, the effect of the inelasticity 
is easily incorporated by changing the binary collision rule. As in the elastic case, the 
dynamics of the system consists of a sequence of free streaming between collisions and 
instantaneous velocity changes on binary collisions. Therefore, a statistical mechanics 
description can be formulated by means of the Liouville equation of the system [21 E|. 
From a conceptual point view, this description has the same degree of justification as 
that of inherently non-equilibrium molecular systems. Proceeding to a more mesoscopic 
description, as provided by kinetic theory, is a more difficult and controversial problem. 
Formally, the same methods as developed for hard sphere gases with elastic collisions 
can be applied. In particular, in the low density limit, the inelastic Boltzmann equation 
is derived. Thus, it seems that the basis for this equation is the same as in the elastic 
case. A different question is the domain in the parameter space of density and coefficient 
of restitution in which it actually provides an accurate description of the system [3] . 

The Boltzmann equation is a kinetic equation describing the time evolution of 
the one-particle distribution function, that is an average quantity. It is also possible 
to studyagain in low density limit, fluctuations and correlations in the positions and 
velocities of the particles by means of the associated distributions functions. It is known 
that correlations can become quite large in out of equilibrium elastic systems, on the 
appropriate long time and length scales. This should not be interpreted as a signature 
of the failure of the Boltzmann equation. The latter only requires that the precollisonal 
part of the correlations of two particles at contact be small E]. Closed kinetic 
equations for the fluctuation and correlation functions valid in the low density limit, can 
be derived inside the same approximation schemes that lead to the Boltzmann equation. 
This includes formal expansions in a small parameter |7j and cluster expansions of the 
reduced distribution functions [8j. Both methods have been considered for granular gases 
[HI EI] ■ Nevertheless, while the inelastic Boltzmann equation and its implications have 
been extensively studied, little is known about fluctuations and correlations in granular 
gases, even though it is widely recognized that they are crucial for the understanding of 
the properties of these systems. 

As mentioned above, granular gases display fluid dynamical behavior. Then, they 
are expected to obey some kind of hydrodynamic equations, understood as a set of 
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closed equations for the hydrodynamics fields: density, flow velocity, and (granular) 
temperature. In fact, more or less phenomenological hydrodynamic equations have 
been successfully applied to many different situations. An appropriate context in which 
the issue of the existence of hydrodynamics for dilute granular gases, the form of the 
hydrodynamic equations, and their range of validity, can be addressed is provided by 
the inelastic Boltzmann equation. By using an extension of the Chapmann-Enskog 
procedure, the Navier-Stokes equations have been derived, and formal expressions for the 
transport coefficients obtained. They are, in principle, valid for arbitrary physical values 
of the coefficient of restitution jTTj. Recently, these expressions have been expressed in 
the form of time integrals of the correlations of some velocity functions computed in 
a given reference state, the so-called homogeneous cooling state [121 EH]- Consistency 
requires that these integrals exist in the long time limit. The above expressions are the 
generalization of the Green-Kubo relations for elastic gases, although with significative 
differences that could have been hardly guessed a priori. 

It must be stressed that the Chapmann-Enskog method and the Green-Kubo 
relations presume the existence of hydrodynamics on the appropriate length and time 
scales and, therefore, they do not add any fundamental support to the validity of a 
hydrodynamic description for granular gases. On the other hand, if Navier-Stokes 
equations exist, the Chapmann-Enskog results provide their explicit form. A deeper 
physical understanding of the structure of the Green-Kubo relations and their validity 
follows from the analysis of the linearized Boltzmann equation, and more specifically of 
the spectrum of the associated operator (THj . It is possible to identify the hydrodynamic 
part of the spectrum, that is the one determining the form of the Green-Kubo relations. 
It is given by d + 2 points that in the elastic limit correspond to the d + 2 collisional 
invariants, where d is the dimension of the system. The key remaining point for the 
existence of hydrodynamics is whether the corresponding modes dominate for long 
times and long wavelengths ITH] . Quite interestingly, the hydrodynamic modes 
also determine the form of the fluctuations of the hydrodynamic quantities in granular 
systems, reflecting a generalization of the fluctuation-dissipation theorems. 

The aim of this paper is to review some recent advances and present some new 
results along the above ideas. Particular emphasis will be put on the relevance of the 
hydrodynamic modes for the study of many fundamental and applied issues in granular 
fluids, as it is the case for molecular, elastic systems. Along the paper, and specially 
in the last Section, some of the many remaining open questions in the context of the 
matter discussed here will be indicated. 

2. The homogeneous cooling state of a dilute granular gas 

The system considered is composed of N smooth inelastic hard spheres (d = 3) or disks 
{d = 2) of mass m and diameter a. Their dynamics consists of free streaming until two 



4 



particles % and j are at contact, when their velocities change according with the rule 

Vi -> v i = b a Vi =Vi — (<x ■ Vij) cr. 

I + a ^ ^ 

Vj -> v'j = b^vj = Vj H — (ct • v^) a. (1) 

where = — Vj is the relative velocity, a is the unit vector pointing from the 
center of particle j to the center of particle % at contact, and a is the coefficient of 
normal restitution. It is defined in the range < a < 1 and measures the inelasticity of 
collisions. Here it will be taken constant. 

^From the above dynamics, the Liouville equation governing the time evolution of 
the iV-particle distribution function of the system can be derived 0IH]- An equivalent 
representation is given by the BBGKY hierarchy for the s-particle distribution functions 
(s — 1, 2, . . .). In the low density limit, a closed equation for the one-particle distribution 
function, f(r,v,t), can be derived by using the same procedures as for elastic hard 
particles [21 E] • This is the inelastic Boltzmann equation, 

l i + v--^f(r,v,t) = J[fJ], (2) 

where J[f\f] is the inelastic Boltzmann collision operator |2~1 131 [TGI IT71 ITH] . 

For an isolated granular fluid, no homogenous steady solution is possible. It 
is trivially seen that the energy of the system decreases monotonically due to the 
inelasticity of collisions. On the other hand, there is a spacially homogeneous state, the 
so-called homogeneous cooling state (HCS), for which all the time dependence occurs 
only through the granular temperature. The latter is defined in the usual way from the 
average kinetic energy. For the one-particle distribution, this implies a scaling of the 
form JH] 

f hC s(v,t)=n H v^ d [T hcs (t)]xhcs(c), v [T hcs (t)]= • (3) 

Here Uh is the number of particles density, ks is the Boltzmann constant, and Th cs {t) 
is the time-dependent temperature of the HCS. Finally, Xhcs(c) is an isotropic function 
of c = v/v [T(t)]. Use of Eq. into Eq. (J2J) provides an integral equation for Xhcs, 

~~2~Q^ ' ( C Xhcs) — Jc\Xhcsi Xhcs], (4) 

where J c is the scaled Boltzmann collision operator, 

Ji^d^SfecJrf^l (5) 

T (c, ci) = J da Q[g ■ a]g ■ a [a" a 6~ 1 (c, c x ) - l] , (6) 

with g = c — C\, G the Heaviside step function, and b~ l an operator replacing the 
velocities c and C\ to its right by the precollisional values obtained by inverting the 



5 



collision rules given in Eq. Moreover, the scaling property implies that the time 
dependence of the temperature is given by the Haff law 

T hcs (t) = T hcs (0) ' 



1 + -(hcs(O)t 



(7) 



The cooling rate of the HCS, (hcs, is given by 



(hcs(t) 



1 



n H v (t) 



dc / dcxg^hcsi^Xhcsici), 



and 



Co 



v [T hcs (t)Y (9) 
l/nH<y d ^ 1 is proportional to the mean free path of the particles. Since 

1 /2 

(hcs oc T hcs , Eq. ((7|) becomes independent of the initial condition in the long time limit, 



where 



T, 



hcs 



4(a)- 



c 



(hcs (t) 



(10) 



indicating that all the homogeneous cooling states of a given system tend to converge in 
the long time limit 30J. This seems to suggest the existence of some kind of H theorem 
for the homogeneous Boltzmann equation for inelastic hard particles, although no proof 
of such a theorem has been given up to date. 



3. Linearized Boltzmann equation and hydrodynamic modes 

Now, the initial value problem for small perturbations around the HCS of an isolated 
granular gas will be considered. In this case, the Boltzmann equation can be linearized 
around the HCS, and it will be shown that, in appropriate dimensionless units, its 
general solution can be related to a linear eigenvalue problem. By considering times 
over which the deviation of the distribution function from fh cs remains small, we can 
write 

f(r,v,t) = f hcs (v,t) + 5f(r,v,t), \8f(r,v,t)\^f hcs (v,t). (11) 

In order to eliminate the time dependence of the HCS, it is convenient to introduce the 
scaled velocity c defined below Eq. and the dimensionless quantities 

I=j, s = j\t' V -^, 5 X (l,c,s)=n- H 1 v d (t)5f(r,v,t). (12) 

The time scale s is proportional to the average of the accumulated number of collisions 
per particle in the interval (0,t). Then, substitution of Eq. (jllj) into the Boltzmann 
equation and retaining only terms up to first order in the deviations, gives the linear 
Boltzmann equation 

d d 

■^~ 5 X + c- t^5x = A5x(l,c,s), (13) 
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where the linear operator A is given by 

Co 9 

= J c [Xhcs, Sx] + M 5 X, Xhcs] ~ y 7^ • (cSx) ■ (14) 

In this representation, the effect of the cooling appears explicitly through the last term 
on the right hand side of the above expression. 

Solutions to Eq. (fT3~J) will be sought in a Hilbert space defined by the scalar product 

<g\h>= J dc X £(c)g*{c)h(c), (15) 

where g* denotes the complex conjugate of g. It must be noted that this is a stronger 
requirement than the necessary normalization of the distribution function. Since the 
collision operator A does not affect the space variable, it is enough to consider a single 
Fourier mode, i.e., 

6 X (l,c,s) = e lk4 6x(k,c,s), (16) 
so that the linearized Boltzmann equation becomes 

^-5 x =(A- t k-c)5 X . (17) 

All linear excitations are determined from the spectrum of the operator A — ik ■ c. 
The hydrodynamic part of this spectrum is defined as follows. For k = it is given 
by those eigenvalues that coincide with the eigenvalues of the balance equations for 
the momentum density, flow velocity, and temperature obtained from the linearized 
homogeneous Boltzmann equation. Assuming that the cooling rate for homogeneous 
systems can be expressed as a function of the density and the temperature, these 
eigenvalues are found to be 

Ai=0, A 2 = Co/2 A 3 = -Co/2, (18) 

the eigenvalue A 2 being <i-fold degenerated. For finite k, the hydrodynamic modes are 
defined as the solutions of the equation 

(A-ik-c)Uk,c) = \&(k,c) (19) 

that are continuously connected as functions of k to the above special solutions for 
k = 0. Of course, consistency requires that A includes the spectrum of the balance 
equations for k = 0. Straightforward calculations show that jT^j 

A&(c) = Ai&(c), i = 1,2,3 (20) 

with the Aj's given by Eq. (j!8)l and 

d 

6(c) = Xhcs(c) + — ■ [cXhcs(c)] , 

<hh ' Jr) &(c) = -#-[c Xto (c)]. (21) 



dc dc 

This confirms the existence of the hydrodynamic modes. Their evaluation for finite k 
is now a technical problem. To Navier-Stokes order, they have been computed by using 
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perturbation expansion methods [14], and the results agree with those obtained from 
the Navier-Stokes equations derived by means of the Chapmann-Enskog procedure jllj . 
Here an alternative approach, designed to derive in a direct way expressions for the 
Navier-Stokes transport coefficients, will be outlined. 

As a consequence of A being not Hermitian, the eigenfunctions £j are not orthogonal. 
It is then useful to introduce a set of functions {£;} being biorthogonal to the set 
i = 1,2,3}. A convenient choice is 

{&} = ^Xhcs(c),c X hcs(c), + Xh cs (c)| , = 5ij . (22) 

The general solution to the linearized Boltzmann equation in the Hilbert space can 
be formally written as 

5x(k, c,s) = J2 a i(k, s) Xi (k, c) + 6x m (k, c, s), (23) 

i 

where the sum extends over the d + 2 hydrodynamic modes and the term 5x m contains 
all the other "microscopic" excitations. Assuming that the hydrodynamic spectrum 
dominates for long times and small gradients, Sx m can be neglected, and the coefficients 
5j are given by 

~a i (k,s) = &\5x(k,s)). (24) 

These coefficients are combinations of the hydrodynamic fields, i.e. they can be 
expressed in terms of the deviations of the macroscopic density n, temperature T and 
fluid velocity u. It is now possible to derive a formal expression for 5x(k, c, s) to first 
order in k, and from it the Navier-Stokes order for the pressure tensor and the heat 
conductivity [13J. The results will be given in the next section. 

The above discussion relies on the assumption that the hydrodynamic point 
spectrum found here is isolated and smaller than all the other parts of the spectrum, so 
that it dominates for long times and long wavelengths. This is a necessary condition for 
the existence of hydrodynamics and has been proven in the elastic case. The proof does 
not extend directly to the inelastic case, although there are some reasons to expect that 
this is the case even for rather strong inelasticity. A detailed discussion of the basis for 
a hydrodynamic description of a granular gas in this context is given in [TH] . 

4. Green-Kubo representation of the transport coefficients 

The procedure described in the previous Section leads to the following expressions for the 
dissipative part of the pressure tensor Hij(r, t) and the heat flux q(r, t) to Navier-Stokes 
order [13]: 

% = -V + ^ - -foV ■ u) , q = -kVT - ,V«, (25) 
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where 77 is the shear viscosity, k the (thermal) heat conductivity, and fi a new transport 
coefficient that vanishes in the elastic limit and that will be referred to as the diffusive 
heat conductivity. They are given by 

(26) 
(27) 



V = n H m£v {T) / ds e~ s ^ 2 (A xy {s)$ 2jXy ) , 
Jo 

POO 

K = n H £v (T) / dse 8 ^ 2 {Z x (s)$ 3 , x ), 
Jo 



fx = mlvl{T) 



poo -1 poo 

/ ds ( e <o/2 _ X ) (E x ( s )^ 3>x ) + - ds (Z x (s)c x 
Jo 1 Jo 



(2f 



In the above expressions, A xy and are the transversal momentum flux and the heat 
flux, respectively, 



EJc) 



c 2 - 



d+2 



while $2 in and $3 x are "modified fluxes" 



ainxh cs (c)) 



9c 



The angular brackets denote "time-correlation functions" in the HCS, 
(g(s)h)= / dcxhcs(c)g(c,s)h{c) 



e sA{c) g(c) 



where the time dependence is given by 
g(c, s) 

with 

A(c) = / dCiXhcsici) I dd-e(g • a)g • $[b c 



[c,c l )-l](l + P) + ^c~ 



(29) 



(30) 



(31) 



(32) 



(33) 



Here, P is an operator that interchanges the velocities c and C\ to its right. Equations 
(|26 |) - (}2*%j) are the generalization for inelastic hard collisions of the dilute limit of the 
well-known Green-Kubo relations for molecular fluids, reducing to them in the elastic 
limit [TJ]. They can be understood as generalized fluctuation-dissipation relations, 
in the sense that they relate the transport coefficients characterizing dissipation in 
inhomogeneous systems with the time decay of fluctuations in the homogeneous reference 
state (the HCS). Equivalent expressions for the transport coefficients are obtained by 
solving the nonlinear Boltzmann equation by the Chapman- Enskog method [TT] IT2*] and 
also by perturbation analysis of the hydrodynamic spectrum of the linear Boltzmann 
equation ^3]. The effect of the inelasticity of collisions manifests itself in several ways. 
Perhaps the most unexpected among them is that the initial dissipative fluxes A xy and 
S z are replaced by the modified fluxes &2,xy and $3 jX in the time correlation functions 
appearing in the expressions of rj and k. This is a direct consequence of the form 
of the hydrodynamic modes. While in a molecular system they are given by linear 
combinations of the Maxwellian times 1, c, and c 2 , in granular gases they are given by 
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Eq. (J3TJ), and it should be realized that Xhcs(c) is not a Maxwellian outside the elastic 
limit a — 1. 

Since the diffusive heat conductivity /i is a peculiarity of granular systems, it seems 
appropriate to comment about its relevence. In general, its contribution to the heat 
flux is much smaller than the one coming from the usual heat conductivity associated 
with the temperature gradient, because k 3> /i for all a. Nevertheless, there are some 
cases where the existence of this contribution manifests itself in a clear way. Consider, 
for instance, an open granular system in presence of gravity. The system is fluidized by 
means of a vibrating bottom plate. In the steady state, the granular temperature profile 
presents a minimum at a certain height, increasing from there on. This temperature 
inversion has been observed in computer simulations ^Hl EH E] and also in experiments 
[22| 123) E3J EHJ. Its theoretical explanation is tied to the existence of the diffusive 
contribution to the heat flux [T9"j . and the value of the transport coefficient fi can be 
directly obtained from the measurement of the heat flux at the temperature minimum 
[2E]. For a low density gas, the results are in good agreement with the theoretical 
predictions discussed in the next Section. 

5. Simulation study of the Green-Kubo relations 

The formal expressions of the transport coefficients given by Eqs. ()26)) -()28) ) are rather 
involved and have been evaluated analytically only in some approximation, namely by 
introducing an expansion in orthogonal functions, keeping only rather arbitrarily the 
lowest orders contributions [TTJ EH EH]- On the other hand, the above Green-Kubo 
relations can be transformed into a form that is suitable for evaluation by means of 
iV-particle simulation techniques. Two main points must be addressed for that. First, 
there is the technical difficulty that the operator A defining the time dependence of 
the fluxes, involves the cooling rate Co that, therefore, should be known a priori and it 
is not. This can be overcome by making a change in the time scale and realizing the 
existence of an exact mapping of the HCS onto a steady state [23 [3Dj- The second 
issue is to transform Eqs. 1)26)1 - (|2"K|) into expressions involving a well-defined iV-particle 
dynamics, instead of the effective Boltzmann dynamics defined by A. This can be done 
by using the same kind of assumptions as needed to derive the Boltzmann equation and, 
in particular, assuming that velocity correlations in the HCS are negligible in the low 
density limit 30J. 

We have used the direct simulation Monte Carlo method (DSMC) 31J to simulate 
the iV-particle dynamics of a dilute granular gas composed of inelastic hard spheres 
(d = 3). The details of the practical implementation of the above ideas and the 
simulation technique are discussed in refs. (SOj and [321 ■ Here it will only be 
stressed that the DSMC allows to simulate a homogeneous state in such a way that 
it stays homogeneous "by definition", avoiding the spontaneous developing of spacial 
inhomogeneities. This is important in the present context, since the HCS is known to 
be unstable against spacial long wave perturbations, even in the low density limit. This 
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Figure 1. Left: The dimensionless reduced coefficient of shear viscosity 77* as a 
function of a. The symbols are from the DSMC method, while the solid line is the 
analytical result obtained in the first Sonine approximation Right: The same for 
the reduced coefficient of thermal heat conductivity k* 



is the so-called clustering instability [SSI Ell- 
in the simulations, it is observed that the HCS, or more precisely its equivalent 
steady representation, is reached after a few collisions per particle. This supports that 
the hydrodynamic point spectrum identified in Sec. E]is isolated and smaller than all the 
other parts of the spectrum. Moreover, all the time correlation functions appearing in 
Eqs. (J2EJ)-(|2EJ) are found to decay exponentially, at least until they decay by two orders 
of magnitud from their initial values. This has been exploited in order to compute the 
time integrals. The results obtained for the shear viscosity are shown in Fig. Q (left). 
What is actually plotted is the shear viscosity scaled with its elastic value in the first 
Sonine approximation, i.e. 



V = 



rj(T) 



5(mk B T) 



1/2 



(34) 



Also plotted in the same figure is the approximate theoretical prediction obtained when 
computing Eq. (|2fi|) in the first Sonine approximation ^Tj. It is observed that the 
viscosity monotonically increases with the inelasticity of collisions (decreasing a) and 
that the theoretical prediction is in good agreement with the simulation results over all 
the range of values of a. 

The simulation results for the (thermal) heat conductivity are given in Fig. ^ 
(right). Similarly to above, the dimensionless plotted quantity is 



K 



k{T) 75k B (k B T)^ 



(35) 



Ko(T)' 64(vrm) 1 /V 2 ' 

Again, it is found that the transport coefficient is a monotonic decreasing function of 
a. However, in this case the first Sonine approximation leads to significant quantitative 
discrepancies for small values of a. This may be, at least partially, due to the fact that 
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V- 1.0 




Figure 2. The same as in Fig. ^ but for the reduced coefficient of diffusive heat 
conductivity (i* 



the fluxes appearing in the expression of the heat conductivity involve higher velocity 
moments than those present in the expression of the shear viscosity. 

The dimensionless reduced diffusive heat conductivity /i* has been defined as 
* _ nfi(T) 

" = tWV (36) 

and the results obtained for it are presented in Fig. |21 As it is the case for the 
thermal heat conductivity, the first Sonine approximation significantly overestimates 
this transport coefficient for strong dissipation. 

The analysis of the transport coefficients for a dilute gas of inelastic hard disks 
(d = 2) has been reported in [32]. The results are quite similar, the main difference 
being that the discrepancies of the first Sonine approximation predictions are stronger 
for d = 2 than for d — 3, as it is usually the case. 

As already pointed out, in order to transform the effective one-particle dynamics, 
implicit in the linearized evolution operator A, into an iV-particle dynamics, it has been 
assumed that velocity correlations are negligible in the HCS of a low density granular 
gas. Nevertheless, when the DSMC data are carefully analyzed, it is seen that for strong 
inelasticity (a < 0.6), velocity correlation effects are clearly identified F° r instance, 
their contribution to the initial value of the time-correlation functions appearing in the 
iV-particle representation of the Green-Kubo expressions can be up to 10% of the total 
value. This raises some fundamental issues. Of course, a first possibility is that the 
velocity correlations be just an artifact introduced by the simulation method itself, i.e. 
an inaccuracy of the DSMC method. On the other hand, if they are not, how can these 
correlations be incorporated in a theoretical description of dilute granular gases? Is it 
consistent with the description provided by the inelastic Boltzmann equation? A partial 
answer to these questions will be given in the next Section. 
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6. Velocity correlations and total energy fluctuations in the HCS 

In spite of the great advance in the understanding of rapid granular flows in the last 
years, little is known about fluctuations and correlations in those systems. Most previous 
studies have focussed on the study of the initial build-up of spacial correlations in the 
development of the clustering instability [HHJ EE] or on the short time and length scales 
precollisional velocity correlations when starting from a thermal equilibrium state j37j . 
Recently, an empirical expression for the two-particle distribution function of a system 
in the HCS, based on simplicity and symmetry considerations, has been proposed [38J. 
It includes both spatial and velocity correlations, but their relation with the dynamics 
of the system has not been analyzed. 

Closed equations for the distribution functions characterizing correlations in a dilute 
granular g be derived by using the same kind of approximations as needed to 

derive the inelastic Boltzmann equation. In particular, in the hierarchy method it 
is assumed that the order in the density of the precollisional part of the correlation 
functions of two particles at contact, increases as the number of involved particles 
increases jHj- Of course, the results obtained under this assumption must be consistent 
with it. 

Here we will restrict ourselves to a dilute gas in the HCS. It is assumed that the 
Liouville equations admits a solution describing this state and having the scaling form 

Pke.{{ri,Vi},t) = [v Q {t)r dN Plc S {{rij},{<k}), (37) 

that is invariant under space translations. This implies an analogous scaling for all 
the distribution functions of the system and, in particular, leads to Eq. (JHj) for the 
one-particle distribution function. Two-particle correlations are described by means of 
the function ghcs( r i2, v i, v 2, t), defined through the cluster expansion of the two-particle 
distribution function f2,hcs(f 12, v%, v 2 , t), 

f2,hcs( r 12, Vi, V 2 , t) = fhcs{Vl,t)f hcs (v2, t) + 9 hcs{r 12, «1, «2, *)• (38) 

As a consequence of the scaling property (|37j). the reduced distribution function 

9hcs(h2, ci, c 2 ) = n~ H l l d vl d {t)g hcs (ri2, v u v 2 , t) (39) 

does not depend on s. Then, from the second equation of the BBGKY hierarchy, 
neglecting the term involving three-particle correlations, a closed equation for g^cs is 
found P33, 



d 

C12-TT, A(ci) - A(c 2 

Oils 



9hc S (h2, Cl, c 2 ) = 5(l 12 )T (c 1 , c 2 )xhcs(c l )xhcs(c 2 ). (40) 



In order to focus on the velocity correlations, we introduce the marginal velocity 
correlation function 

iphcs(ci, c 2 ) = J dl 12 g hcs (l 12 , ci\ c 2 ) +5(c 1 - c 2 )xhcs{ci)- (41) 

Integration of Eq. (jjQj) yields 

[A(ci) + A(c 2 )] Vao»(ci, c 2 ) = -T (ci, c 2 )xhcs(ci)Xhcs(c 2 ). (42) 
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Trying to find the complete solution of the above equation is a formidable task. We 
will restrict ourselves to find the hydrodynamic part ip h h } s of the solution, defined by 

^£(ci, c 2 ) = Viphcs(c h c 2 ), (43) 

where V is the projection operator onto the direct product of the hydrodynamic part of 
Hilbert space. Next, we apply V to both sides of Eq. (jUJ) and, in order to get a closed 
equation for ip^ s , we make the approximation PA(cj) = VA(ci)V. Then, the resulting 
equation can be solved yielding [TU] 

VS(Cl,e2) = 083(0)6(^)6,(02). (44) 

The coefficient 033(a) is a functional of Xhcsi vanishing in the limit a — > 1. An explicit 
expression can be obtained in the first Sonine approximation jTU] . The fact that there are 
no contributions to i/j h h J s from the other combinations of the hydrodynamic modes is a 
consequence of the symmetry of the HCS and the strict conservation of the number 
of particles and the total momentum. On the other hand, Eq. (|44jl clearly shows 
the presence of velocity correlations that are due to the inelasticity of collision being, 
therefore, inherent to granular systems. 

A direct test of the accuracy of the above expression by means of A^-particle 
simulation is beyond the present possibilities. Nevertheless, it is possible to carry out 
an indirect check by computing the fluctuations of the total energy E. They are given 
by 

([5E(t)}\ cs = (E 2 (t)) hcs - (E(t))l s 

= Nk 2 B T 2 hcs (t) J dc x J dc 2 clcl^( Cl ,c 2 ) 

= Nk 2 B T 2 cs (t)e(a), (45) 

with e(a) = d 2 a 33 (a). In the above expressions, (■ ■ -)h C s denotes ensemble average with 
Phcs({ri,Vi},t). 

To put the above expression in a proper context, it is important to realize that, 
although its derivation relies on the knowledge of the hydrodynamic modes, the validity 
of the result is not restricted to any kind of hydrodynamic regime. On the contrary, it is 
a general result for the HCS, assuming it exits, with the only operational approximation 
given above Eq. (|44jl . We have measured the energy fluctuations by performing 
Molecular Dynamics simulations of a freely evolving system of inelastic hard disks. 
Periodic boundary conditions were used, and the system size was always chosen smaller 
than the critical value needed for the spontaneous formation of spacial instabilities [33J. 
The results obtained are given in Fig. |3]for two small densities, namely uh = 5 x 10~ 3 o~~ 2 
and tih = O.Olcr -2 . The plotted quantity is a\ st = ((5E) 2 )hcs/Nk B T 2 cs) which does 
not depend on time, and each of the reported values corresponds to an average over 
several trajectories. Also plotted is the theoretical prediction provided by Eq. (}45|) . 
i.e. the function e(a) computed in the first Sonine approximation. A quite good 
agreement is observed over all the range of a considered, specially at the lowest density. 
This was expected given the low density character of the theoretical approach sketched 
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Figure 3. Scaled second moment of the energy fluctuations a% st as a function of 
the restitution coefficient a. The symbols are from MD simulations of a system 
of hard disks and the solid line the theoretical prediction given by Eq. I|45|). Left: 
n fl = 5x lCr 3 <j- 2 ; right: n H = O.Olcr -2 

above. Unfortunately, considering smaller values of the restitution coefficient requires 
the simulations of too small systems in order to avoid the instability mentioned above, 
and consequently finite size effects may become relevant. 

Then, we conclude that there are intrinsic velocity correlations in a granular fluid as 
a consequence of the inelasticity of collisions, even in the low density limit. Moreover, 
they increase as the value of a decreases, as expected. Also, the results in Fig. |3] 
indicate that what we have termed hydrodynamic part of the velocity correlations is 
well described by a consistent low density theory, at least for a > 0.7. On the other 
hand, it is not possible to infer anything more precise from the results in this Section 
with regard to the manifestation of velocity correlation effects in the evaluation of the 
Green-Kubo expressions, as discussed in the previous Section. There are two reasons 
for this: 1) it is easily realized that the velocity correlations involved in the Green-Kubo 
expressions lie outside the hydrodynamic subspace, and 2) as already indicated, their 
effect become appreciable for a < 0.6, i.e. for systems more inelastic than those consider 
in the MD simulations reported here. 

^From the MD simulation results, it is also possible to compute the normalized 
distribution of the scaled energy fluctuations SE/(E)h cs ■ An example, corresponding 
to a = 0.8 and uh = 0.02a -2 , is shown in Fig. 0J The results are well fitted by 
a Gaussian, within the numerical precision of the data. In particular, no asymmetry 
is observed. Up to date, there is no theoretical explanation available for this simple 
behavior in an inherent non-equilibrium system. 

7. Discussion 

In this paper, several questions related with the description of dilute granular gases 
have been addressed. The basic theoretical ingredient for the discussion has been the 
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Figure 4. Normalized distribution of the relative energy fluctuations for a system of 
inelastic hard disks with ri# = 0.02er~ 2 . The symbols are from MD simulations and 
the solid line is a fit to a Gaussian. 

identification of the hydrodynamic part of the spectrum of the linearized Boltzmann 
equation for smooth inelastic hard spheres or disks. In the limit of vanishing wavevector, 
the hydrodynamic modes are determined in a unique way by the exact balance equations 
for the density, momentum, and energy, for all values of the coefficient of restitution. 
Beyond this limit, one has to resort to some approximations, for instance perturbation 
expansions, assuming analyticity and convergence. Moreover, if the hydrodynamic 
modes dominate for long times and long wavelenghts, the hydrodynamic Navier-Stokes 
equations for a dilute granular gas follow, as discussed in Sec. El Therefore, the 
fundamental and controversial question of the existence of a hydrodynamic description 
for granular gases described by the Boltzmann equation can be sharply formulated as: 
are the hydrodynamic eigenvalues smaller than all the other parts of the spectrum? For 
elastic systems, this question has been positively answered jHHl, and it can be expected 
that much more be learned about the spectrum of the inelastic linear Boltzmann 
equation in the near future. 

The procedure indicated above leads to formal expressions for the transport 
coefficients of a dilute granular gas that can be expressed in the form of Green- 
Kubo relations. The energy dissipation in collisions introduces significant differences as 
compared with the usual relations for molecular, elastic gases. The expressions have been 
evaluated by means of an iV-particle simulation technique. The results are in quite good 
agreement with the analytical expressions obtained in the first Sonine approximation for 
not too strongly inelastic systems, but significant discrepancies are observed for small 
values of the coefficient of restitution. An analysis of the simulation data show that they 
are, at least in part, due to the presence of relevant velocity correlations in the HCS, 
even in the low density limit. At what extension this implies that the effect of velocity 
correlations should be incorporated somehow in the theoretical analysis of the transport 
coefficients of a dilute granular gas, is an open question. A related question is the role 
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of velocity correlations for finite density. It has been found that velocity correlations 
might become important quite before a granular system develops significant short range 
spacial correlations [I] . Therefore, extensions of the Boltzmann equation to finite density 
should include the effect of velocity correlations in the precollisonal sphere. This offers 
another challenge to the researchers on granular fluids. 

As a necessary first step in the analysis of the velocity correlations in the HCS 
of granular g distribution function describing partially their presence has been 

computed. Again, advantage has been taken of the knowledge of the hydrodynamic 
modes for vanishing wavevector. The identified velocity correlations have an intrinsic 
inelastic character, vanishing in the elastic limit, although the theoretical analysis carried 
out predicts that they stay small for all values of the restitution coefficient. The presence 
of velocity correlations implies that the total energy of the system fluctuates. We have 
compared the theoretical predictions for these fluctuations with molecular dynamics 
results and found a good agreement. Nevertheless, due to limitations associated with the 
instability of the HCS, the simulation results are restricted to not very strong inelasticity. 
The issue of whether the theory also works in the very small a region deserves future 
work. 
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